Intensity correlations and mesoscopic fluctuations of diffusing photons in cold atoms 
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We study the angular correlation function of speckle patterns that result from multiple scattering 
of photons by cold atomic clouds. We show that this correlation function becomes larger than the 
value given by Rayleigh law for classical scatterers. These large intensity fluctuations constitute a 
new mesoscopic interference effect specific to atom-photon interactions, that could not be observed in 
other systems such as weakly disordered metals. We provide a complete description of this behavior 
and expressions that allow for a quantitative comparison with experiments. 
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A wave propagating in a random medium undergoes 
multiple scattering and the intensity pattern resulting 
from interferences of the scattered waves with each other 
is known as a speckle pattern. The angular and time 
dependent properties of these patterns have been exten- 
sively studied [1-3]. They exhibit coherent mesoscopic 
effects and provide a sensitive probe to scattering proper- 
ties of diffusive media. Quasi-resonant elastic scattering 
of photons by cold atomic gases represents in this context 
an important issue, since it provides a new tool to study 
properties of cold atomic gases such as atomic dynamics. 
Photon propagation in atomic gases differs from the case 
of electrons in disordered metals [4] or of electromagnetic 
waves in suspensions of classical scatterers, due to the ex- 
istence of atomic internal degrees of freedom coupled to 
the photon polarization. Some effects of a Zeeman de- 
generacy on the coherent backscattering [2, 5] have been 
recently investigated in the weak scattering limit [6] in 
terms of phase coherence times [7]. 

The purpose of this work is to study the static angu- 
lar correlation function of photons performing coherent 
multiple scattering in a cold atomic gas. The photon in- 
tensity correlation function between angular scattering 
channels is defined using the transmission coefficient T a \, 
by 
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Here, denotes a configuration average over both the 
position of atoms and their internal degrees of freedom 
(see below) and 8T a b = Tab — Tab- For classical scatterers, 
intensity fluctuations obey the Rayleigh law C a bab = 1. 
In the presence of a Zeeman degeneracy, angular corre- 
lations of speckle patterns and intensity fluctuations be- 
come larger than one. This is a new and genuine meso- 
scopic effect specific to multiple scattering of photons by 
atoms and directly related to interference between am- 
plitudes associated to different atomic quantum states. 

Atoms are modeled as degenerate two-level systems de- 
noted by \jgTrig) for the ground state and \j e m e ) for the 
excited state, where j is the total angular momentum and 



m is its projection on a quantization axis. The levels are 
degenerate with \m g \ < j g and |m e | < j e . 

We refer to the following possible experimental setup 
(Fig.l). A light pulse is incident along a direction s a 
onto a dense enough atomic gas confined in a slab ge- 
ometry. This pulse is detected along a direction Sj, after 
being multiply scattered (ab channel). A time r later, 
a second pulse that corresponds to the a'b' channel, is 
detected. We assume that the time r is short enough so 
that the atoms stay at rest between the two pulses. The 
same measurement is repeated after a time r>T, dur- 
ing which the scatterers move. The averaging over spa- 
tial disorder results from this motion. This is a Young- 
like experiment. Thus, although a pulse contains many 
photons, the transmitted intensity T a \, is proportional to 
the probability of one "representative" photon incoming 
along §„, to emerge along s b . 




FIG. 1: Photons in each pair of pulses are scattered by 
atoms at identical positions n but with distinct and uncorre- 
cted quantum numbers (rrn,mf) and (jn'i,m'f). After a time 
T ^> t, a new measurement is performed after the atoms 
moved to new positions r£ . 

The average transmission coefficient T a b is obtained by 
summing all the possible scattering amplitudes, An R ' m ^ , 
corresponding to a given configuration {R, m}. Here {R} 
accounts for the spatial positions of all scatterers, and 
{m} is a notation for their internal Zeeman states both 
before and after scattering. The index n denotes one 
possible multiple scattering path. Squaring the sum of 
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amplitudes we have [8] 
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where TTT denotes a configuration average, over both {R} 
and {to}. When averaging over {i?}, all cross terms 
n / n' vanish because of large fluctuating phase shifts, so 

thatr ab = £J^ m} | 2 - This expression, known as the 
intensity Diffuson, is the leading approximation in the 
weak disorder limit k$l 3> 1, where fco and I are respec- 
tively the wave number and the elastic mean free path 
of photons. The (ab) pulse contains many photons, and 
each of them may change the internal state of atoms. 
Therefore, if {to} (resp. {to'}) is the atomic internal 
configuration seen by a "representative" photon of the 
ab (resp. a'b') pulse, then we can assume that there is no 
correlation between {to} and {to'}. 

Similarly to the average intensity (2), the correlation 
of the transmission coefficient is 
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since, as before, the averaging over {R} leaves only pairs 
of amplitudes having exactly opposite phase shifts. To 
leading order in the weak disorder limit, the only non 
vanishing contributions involve two Diffusons, i.e. two 
possible pairings of amplitudes, cither i = j, k = I, which 
gives T a bTa'b', or i = l,j = k so that 
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The correlation function thus appears as products of two 
amplitudes, that correspond to different internal con- 
figurations {to} and {to'}, but to the same scattering 
path i (or j). Most of multiple scattering paths i and 
j do not share common scatterers so that we can aver- 
age A\ m ^ A\ m and A^ m ^ A^ m ^* separately, since these 
averages are taken upon different atoms, and finally, 



(5) 



In the theory of multiple scattering it is helpful to use a 
continuous description [2] . In this framework, one defines 
two Diffuson functions V^ 1 ^ by [8] 
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and 
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where As a: h = s a , b — s a > tb i. The intensity Diffuson 

p( l )(r, r') is the sum of the terms |A^ m ^(r, r')| 2 be- 
tween endpoints r and r'. On the other hand, the 
correlation Diffuson (r , r' ') is the sum of the terms 

A\ m \v, v')A\ m **(r, r'), i.e., that involve uncorrelated 
configurations {to} and {to'}. 

The two functions 2?( 1 > C ) are obtained from the itera- 
tion of a proper elementary vertex y( l ' c \ that describes 
the microscopic details of the scattering process. The it- 
eration of the elementary vertex is written symbolically 
(either for £>(v) y(*,c) we ghall denote by V y} as 



V = V + VWV + 
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The term V accounts for a single scattering and VWV 
represents its iteration. The quantity W describes the 
propagation of the photon intensity between successive 
scattering events and it will be described later on. 

The elementary vertex is obtained by the pairing of 
two scattering amplitudes of a photon by an atom. It is 
given by [9] 

V = 2^+1 ^^B m 2\V(ii,i2)\jgmi)(j g m 4 \V(i3,i4)\jgm 3 )* 
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where the operator V(i', i) = J2 m ' ( ^\je'm e )(j e m e \d-i 
results from the dipolar interaction energy — d.E between 
atoms and photons, d being the atomic dipole operator 
and E the electric field operator. In the case of V^°\ 
each one of the two coupled scattering amplitudes in (9) 
might belong to a distinct atomic configuration, mean- 
ing that we must consider two distinct couples of initial 
{\jg m i)i \jg m 3)) an d final (\jgin2), \jgTn4}) atomic states, 
as well as two initial (£1,63) and final (£2, £4) photon 
polarization states. The summations over the quantum 
numbers mj result from averaging over initial atomic 
states and from non detected final states. Thus cor- 
responds to the most general case regarding the rrij quan- 
tum numbers. In contrast, corresponds to the differ- 
ential scattering cross-section, for which we set to-i = TO3, 
to 2 = TO4, i\ = £3 and i 2 = £4 in (9). This is because the 
intensity Diffuson is built out of two coupled amplitudes 
that must belong to the same scattering process. This 
distinction between and V^°\ and therefore between 
T>^> and T>^ , occurs only for j g > and it is at the basis 
of the new results we obtain here for mesoscopic speckle 
correlations. 

The iteration (8) is implemented using the decomposi- 
tion of the various terms into standard basis components, 
thus leading to the definition of rank four tensors such as 
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Likewise, the iteration equation (8) for the Diffusons ac- 
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quires a tensorial structure, which reads 



identify the set of characteristic times 
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Here VV = Wb, the function VF describes the scalar part 
of the photon intensity propagator and b, defined by 

b a f3, 7 s = ((S ai - (-) 7 s Q s_ 7 ) (S 05 - (-) /3 s_ /3 s (5 )}, (12) 

accounts for the polarization dependent part. This ex- 
pression follows at once by noticing that after being scat- 
tered by an atom, the two outgoing photon amplitudes 
propagate with a wavevector s = k/fco, random in direc- 
tion but identical for both, and with two different po- 
larization components. Since s is random, the intensity 
propagation is averaged (• • •) over photon wavevectors 
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For (mi = rri3, mi = m.4 in (13)), it is straightfor- 
ward to check that (V^b) a /3^s admits a spectral decom- 
position over 3 projectors only [2, 6, 12]. In con- 
trast, for there are no constraint on the quantum 
numbers, so that the total angular momentum needs not 
to be conserved, and the corresponding spectral decom- 
position involves usually more than 3 projectors T^ K \ 
For a non degenerate ground state level (j g = 0), angu- 
lar momentum is automatically conserved and the two 
vertices become identical, = V^. 

The poles that occur in (15) correspond to diffusive 



direction. The term 5^ - {-) v s^S- v expresses transver- modes of lifetime T £.<0_ This shows up when rew riting 



sality. 

To proceed further, we use the Wigner-Eckart theorem 
to rewrite the tensor V a /3 n s in terms of a summation of 
product of 3j-symbols, 
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(14), with the help of (10), in real space 
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where 



where aj g j e — (2j e + l)/3(2j g + 1). The two tensors and a corresponding expression for Y$ obtained by set- 
ting a' = a, b' = b and in the previous relation. 
The scalar Diffuson propagator V(r,r',t) obeys a diffu- 
sion equation whose solution for a slab of width L is well 
known [2, 10] and leads for (6) and (7) to 



fra/3,7<5 and V Q/ 3. 7< 5 can be written in the form of a 9 x 9 
matrix. According to the spectral decomposition theo- 
rem, they can be decomposed using an orthonormal set 
of (generally) nine projectors T^ K ^ [10]. Looking at (11), 
we wish to find the spectral decomposition of V using 
the spectral decomposition of V and bV. The problem is 
that they do not share the same projectors set in their 
spectral decomposition. We are thus led to define a new 
tensor U by V = UV. It obeys the iteration equation 
U = 1 + WUVb and it involves only the spectral de- 



composition of V^'^b = J2 K=0 
immediately to 
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This leads 
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with V^ c) = T^V {l ' c) and 
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where q (with q = |q|) is the Fourier variable of the differ- 
ence R = r' r between the two endpoints of a multiple 
scattering sequence. The r.h.s in (15) is obtained by us- 
ing the diffusion approximation (i.e. ql <C 1), so that 
W(q) ~ 5- ' 



jr— (1 - q I /3), where D = ci/3 is the pho- 
ton diffusion coefficient [2] and c the speed of light. We 



^=£%>(^(0)) (19) 
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and 
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We have defined the quantities q p = koAs aj Qk' c \x) — 
^Jx 2 + (1/Dt£' c) ) and V{x) = smh 2 (xl)/(xlsmh(xL)). 

We now analyze expressions (19) and (20), which con- 
stitute the main results of this paper. First, consider 
the modes of the average intensity Z>W. It is easy to 
check that is infinite as a result of the Ward iden- 
tity — fa.,-^ (see (16)). The corresponding Gold- 
stone mode Uq^ oc 1/Dq 2 expresses energy conservation 
and long-range propagation of the average intensity. The 

(i) 

two other modes U K are exponentially damped with the 

times Tj£ (see (17)). This expresses photon depolariza- 
tion in multiple scattering [2, 6, 7, 11]. 
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The spectral decomposition of T>^ gives rise to nine 

(c) 

modes and their corresponding times t k ' . Such times 
are well-known to occur in quantum mesoscopic physics 
e.g. in conductance fluctuations of disordered metals in 
the presence of magnetic impurities [2, 13]. The times 
t$ describe how underlying interferences between mul- 
tiply scattered waves (electrons in metals, photons in the 
present case) are washed out in the presence of other 
degrees of freedom. The surprising and new feature of 
the atom-photon scattering, is the occurrence of a mode 
(K = 0) with a negative Tq C L According to (17), this cor- 
responds to an amplified mode that enhances the angular 
correlation function. This amplified mode is present for 
a degenerate atomic transition (j g ,j e > 0) and vanishes 
otherwise. Its origin can be traced out from the vertex 
(9) which can be written as a sum of an incoherent con- 
tribution (i.e., a sum of probabilities) present both in 
V*- 4 -* and V^ c \ and a coherent contribution (i.e., a sum 
of products of quantum amplitudes associated to inter- 
ferences between different Zeeman states). The coherent 
contribution enhances with respect to and its 
iteration gives rise to the amplified mode characterized 
by Tq C \ Expressions (19) and (20) lead to an expres- 




FIG. 2: Angular correlation function C a t a 'b' plotted as a func- 
tion of q = fc()As a for an atomic transition between non de- 
generate (j g — 0,j e = 1, dashed line) and degenerate energy 
levels (j g — l,j e — 2, solid line). The non degenerate case 
describes the clasical Rayleigh scattering of a polarized wave. 
These curves correspond to L — 7l e . The inset gives the de- 
pendence of intensity fluctuations C a bab upon the width L (the 
dashed line is the Rayleigh law). 

sion of C a ba'b' plotted in Fig. 2. This expression allows to 
recover the limiting case of a scalar classical wave [1] cor- 
responding to Y~( c ) = yW = l and to a single mode with 
infinite Tq' c ^ . It also provides a simple expression for an- 
gular speckle correlations of classical Rayleigh scattering 
of a polarized wave [10, 14]. For atomic transitions be- 
tween degenerate levels, Caba'b 1 exhibits a steep decrease 
and large intensity fluctuations (see Fig. 2) as compared 
to the case of non degenerate atomic levels. Intensity 
fluctuations measured by C a bab become larger than one, 
unlike Rayleigh law, C a bab — 1, well-obeyed by classi- 
cal scatterers. Large intensity fluctuations result from 

(c) 

the amplified mode Tq which leads to a divergence of 
the integral in (17). This divergence is cutoff by other 



dephasing mechanisms, such as Doppler shift, inelastic 
scattering or finite size of the atomic trap. Denoting by 
A this upper cutoff, and as sumi ng that the dominant con- 
tributions to T a b and to ST 2 , are given respectively by 
the Goldstone and the amplified modes, we deduce from 
(17) that 

TfT ( xad/l 2 _ i\ 2 
Cabab = -^^A^ (- -J (21) 

where X = (L/L^) 2 — ir 2 involves the diffusion length 

= \J D\t^\ and A is a constant equal to one for 
a non degenerate ground state. This approximate ex- 
pression reproduces the main features of C a b a b plotted 
in the inset of Fig. 2, namely that it can be larger than 
the Rayleigh term and that it is peaked at a value of L 
that depends on the cutoff A. When and A arc in- 
finite, C a bab becomes independent of L and is given by 
the Rayleigh law. Relative fluctuations as given by (21) 
thus provide a direct probe of dephasing mechanisms in 
cold atomic gases. 
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